Impure Thoughts on Inelastic Dark Matter 
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The inelastic dark matter scenario was proposed to reconcile the DAMA annual modulation with 
null results from other experiments. In this scenario, WIMPs scatter into an excited state, split from 
the ground state by an energy 5 comparable to the available kinetic energy of a Galactic WIMP. We 
note that for large splittings 5, the dominant scattering at DAMA can occur off of thallium nuclei, 
with A~205, which are present as a dopant at the 10"'^ level in Nal(Tl) crystals. For a WIMP mass 
m^ « 100 GeV/c'^ and 5 ~ 200 keV, we find a region in 5 — tti^,.— parameter space which is consistent 
with all experiments. These parameters in particular can be probed in experiments with thallium 
in their targets, such as KIMS, but are inaccessible to lighter target experiments. Depending on 
the tail of the WIMP velocity distribution, a highly modulated signal may or may not appear at 
CRESST-II. 



Introduction: For more than a decade, the DAMA 
collaboration has employed ultra-pure Nal(Tl) crystals to 
search for dark matter scattering off a laboratory target. 
Their observation of a modulation in the spectral rate at 
low energies [1] is a challenge to understand. While no 
obvious background can mimic this modulation, conven- 
tional models of dark matter explaining this modulation 
predict signals which would have long since been seen at 
other direct detection experiments [2-8] . In light of these 
tensions, various proposals have been put forward to ex- 
plain the DAMA modulation, such as light dark matter 
with or without ion channeling [9-14], spin-dependent 
scattering [15-18], mirror dark matter [19], momentum- 
dependent scattering [20, 21] and inelastic dark matter 
(iDM) [22, 23]. 

In the iDM framework, WIMPs (weakly interacting 
massive particles) with mass m^^ scatter only by tran- 
sitioning to a heavier WIMP state x* ^ with mass split- 



ting 5 = m^ 



liv comparable to the available 



kinetic energy, which depends on the reduced mass fi. 
Among other things, this kinematical change pushes the 
expected signal to higher energy, increases the modula- 
tion amplitude, and favors heavier target materials. To- 
gether, these features allow a positive signal at DAMA 
while suppressing or eliminating signals at other experi- 
ments [24-27]. However, recent results have placed iDM 
under increasing pressure. The null results at CRESST- 
II [5, 28], CDMS [3] and ZEPLIN-III [29] limit the al- 
lowed parameter space to non-Maxwellian halos [30, 31]. 
While such properties may be natural, it is clear that 
iDM is now very constrained. As such, it is worth revis- 
iting the original proposal to consider whether important 
effects may have been neglected. 

Inelastic Sensitivities: Owing to the introduction of 
the splitting parameter 5, the kinematical requirement 
for scattering becomes 



where tun is the mass of the target nucleus, fj, the re- 
duced mass of the WIMP-nucleus system and En the re- 
coil energy. Due to this constraint, different target nuclei 
sample significantly different parts of the WIMP velocity 
distribution. Since the Maxwellian velocity distribution 
is falling exponentially at its tail, different targets with 



different threshold velocities ' 
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can have dramatically 
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different sensitivities to a WIMP with given mass m^ 
and splitting S. As a consequence, even sub-dominant 
components of the target can be the dominant source of 
scattering. In particular, for large values of S, it may be 
impossible for iodine scatterings to occur at DAMA at all. 
However, while iodine is the heaviest element present at 
DAMA in large quantities, the thallium dopant is present 
at the IQ-^ level in the Nal(Tl) scintillator [32]. 

Thallium as a Target: Let us quantify whether the 
dominant scattering in DAMA may be arising from scat- 
tering off of the thallium which is present in the crys- 
tal. The tightest constraint on this high-(5 region of iDM 
parameter space comes from CRESST-II. The atomic 
mass of tungsten (A ^ 184) is nearly as large as that 
of thallium {A ^ 205). Given the small concentration 
of thallium in DAMA, the observed modulation rate of 
^ 0.04 counts/day/kg of Nal(Tl) would correspond to an 
unmodulated event rate of ~ 40 counts/day/kg if thal- 
lium were the target nucleus. As a consequence, the con- 
straints from CRESST-II are that thallium can only be 
the target if scattering off of tungsten is nearly kinemat- 
ically forbidden [49]. 

We show the allowed iDM parameter space from 
DAMA including scattering off of thallium, together with 
the constraints from the tungsten target of CRESST-II in 
figure 1. To this end, we calculate the annual modulation 
at each point in (5 — m^— parameter space and compare 
it to the observed modulation in the binned 2 — 8 keVee 
region from the recent DAMA data [1]. 

We use the Helm form factor [33], with parameteri- 
zation as described in [12], and a Maxwellian velocity 





FIG. 1: The DAMA-allowed range of 5 — tti^— parameter space at 90% confidence for x ~ Tl scattering only (outer green 
hatched region) and constraints from CRESST-II (inner red hatched region). In the remaining allowed range of splittings 5, no 
scattering on sodium or iodine occurs. These contributions to signal at lower S axe neglected here. Left: shown for a Maxwellian 
halo with Hosc = 500km/s, right: using the Via Lactea II simulation data (see [26]). 



distribution with vq = 220km/s and ■^csc = SOOkm/s. 
Wc note that our results do not qualitatively depend on 
this choice: for example, using Wosc = 600kni/s, the re- 
gion looks similar but is shifted by about 40 keV towards 
higher S. We further test the sensitivities by exploring 
other halo models. In particular, we use the tabulated 
Via Lactea II data from [26] , which leads to qualitatively 
similar results, as shown in figure 1. We find the best 
X^ fit parameter point for DAMA in cross section a„, 
splitting S, and mass m^, and consider the 90% confi- 
dence region (or Ax^ < 6.25) for these three parame- 
ters. This region, projected onto the (S, m^) plane, is 
given by the green contour in figure 1. For a given pair 
of {S,m^), we choose the lowest cross section i7„ consis- 
tent with the DAMA data at 90% confidence, and evalu- 
ate the expected CRESST-II signal given this cross sec- 
tion. The cross sections per nucleon near the top of the 
green contour are ^ 2 x lO^'^^ cm^ and higher. While 
these large cross sections are difficult (if not impossible) 
to explain with standard-model mediators, they can be 
achieved with light mediators [34-38]. To place limits, 
shown in red in figure 1, we use the CRESST-II commis- 
sioning run release [5] , and require the signal be less than 
the 7 events observed. As one moves upward in S in the 
allowed range of parameters, while keeping the DAMA 
modulation fixed, the event rate at CRESST-II is drop- 
ping rapidly. Consequently, one could employ more com- 
plicated techniques, such as the maximum gap technique, 
but these would all achieve essentially the same result, 
since over much of the allowed range of parameter space 
to explain DAMA, the scattering signal at CRESST is 
considerably less than 7 events. 

Discussion: There are no measurements of the 
quenching factor for thallium, and thus we must make 
use of estimates. For the DAMA-allowed region, we as- 



sume the quenching q is approximately proportional to 
the path length of the thallium nucleus in the crystal [39] . 
We calculate this path length using the SRIM code [40] 
and find that it scales approximately as mjj . Thus, we 
conservatively take 0.88 > qri/qi > 127/205, where the 
upper limit comes from the ratio of path-lengths using 
the SRIM database, and the lower comes from assuming 
a simple inverse proportionality to mass. Taking a range 
of 0.06 < qi < 0.09 [41-43], we find a range of quenching 
factors 0.037 ^ qn ^ 0.08, and conservatively combine 
all allowed regions in figure 1. Given the uncertainties, 
we opt to include a broader range of quenching factors. 
Quantitatively, our fits prefer lower quenching factors 
(9ti ~ 0.05), but this is a bit misleading. The dominant 
factor contributing to the fit is the location of the form 
factor zero of thallium. For q^i ~ 0.07, this happens to 
fall precisely in the middle of the DAMA energy range, 
using our Helm parameterization. However, this form 
factor has never been measured, and if its zero is shifted 
up even by 5% in momentum transfer q = \/2toti-E_r, 
we can find good fits even with gxi = 0.08. Given these 
uncertainties, we emphasize that only the specifically al- 
lowed range of parameters depends on the quenching fac- 
tor, but not the presence of an allowed range itself. 

The CaW04 crystals employed by CRESST-II contain 
impurities [44], for instance, Nd (~ 1000 ppb), Gd (~ 
4000 ppb) and Sm/Dy/Hf/Os (~ 20 ppb each), but in 
addition to their small abundances, these impurities are 
all too light to be relevant as a target. Heavier impu- 
rities are also irrelevant due to their small abundances. 
Even taking a high rate of ~ 600 ^Bq from ^^°Pb [45], 
the concentration of ^^°Pb compared to "'^^W is ~ 10"^^. 
Concentrations of ^^^U, ^^^U and 232xh [46] can be es- 
timated to be below 10~^^ compared to "'^^W, and thus 



are safely sub-dominant. 

The KIMS experiment uses doped CsI(Tl) targets [6], 
also with a thallium concentration at the 10~^ level [47]. 
However, the event rate ~ 0.28±0.18 counts/day/kg [24] 
observed in KIMS is consistent with a non-zero value. 
Since CRESST-II pushes us to large values of 6, the 
modulation should be 0(1) [24], and the DAM A mod- 
ulation rate of 0.04 counts/day/kg should be interpreted 
as roughly the average annual rate, as well. Thus, KIMS 
does not place a strong constraint. However, it is im- 
portant to note that KIMS should necessarily see the 
modulation, given the doping concentration is similar to 
that of DAMA. 

We have performed the analysis using a Maxwellian ap- 
proximation, but this is well-known to be a poor approx- 
imation for iDM [26] . Consequently, we must be cautious 
in our predictions for CRESST-II. Clearly, the signal at 
CRESST-II must be highly suppressed to compensate for 
the small concentration of thallium at DAMA. Within 
the Maxwellian approximation, it appears that the event 
rate is naturally at higher (~ 40keV) energies. Nonethe- 
less, the presence of high-velocity structures could be 
relevant and cause a non-negligible signal to appear at 
CRESST-II, even at lower energies, although the energy 
range is difficult to specify. Using equation 1, and taking 
a particle with velocity ~ 800km/s, corresponding to a 
particle near the escape velocity after boosting into the 
Earth's frame, we can find the relevant energy ranges al- 
lowed from the parameter space in figure 1. We find that 
recoil energies as low as ~ 15 keV are in principle possible 
at CRESST-II. However, a signal from such high veloci- 
ties would be highly modulated [30] , which would then be 
potentially clear in an annual cycle of CRESST-II data. 

Since the allowed region of parameter space is tuned to 
suppress the CRESST-II signal while leaving a signal at 
DAMA, it should be understood that the precise location 
of the allowed range is sensitive to small changes in form 
factors, halo models, and other inputs, such as the precise 
value of the DAMA doping. That said, the existence of 
a consistent region of parameter space is robust, even if 
its precise position is not. 

Conclusions: We have re-examined the iDM scenario 
keeping in mind the (intentional) thallium impurities of 
the Nal(Tl) crystals used at DAMA. We find this opens 
up a region of parameter space at large splittings S. Con- 
straints from the CRESST-II experiment are strong, and 
the value of S in the remaining region of parameter space 
is tuned at roughly the 5% level. Nonetheless, this region 
exists and is a challenge to probe. A highly-modulated 
signal could arise in CRESST-II, although this is sen- 
sitively dependent on the halo model, and difficult to 
quantify. It is entirely possible that all conventional tar- 
gets are too light to detect such a scenario, making doped 
scintillators the target of choice to probe this region of 
parameter space. In particular, the KIMS experiment, 
with a doped CsI(Tl) target, should necessarily see the 



modulation. Variations in the thallium dopant within the 
DAMA crystals, if controlled, could also test this scenario 
in future expansions of the DAMA experiment. 
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If the dark matter scatters inelastically with spin- 
dependent couplings [27], this CRESST-II constraint 
goes away, since tungsten has almost no spin coupling, 
while thallium has a net proton spin. These scenarios 
would only have signals off of thallium doped scintilla- 
tors. However, the cross sections required by the DAMA 
modulation are quite large, ~ 10"^'^ cm^, which are dif- 
ficult to realize in particle physics models. Inelastic mag- 
netic dipole scattering would be another possibility [48], 
as thallium also has a large magnetic dipole, but again, 
the magnetic moment would have to be roughly larger 
than that of the proton. 



